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I. INTRODUCTION 

Within the circuit model of quantum computing, the 
fundamental building block of a quantum computer is 
the quantum bit or qubit. Our ability to precisely exe- 
cute unitary operations within qubits is an essential pre- 
requisite for the implementation of quantum algorithms 
and for harnessing the full computational power of more 
complex quantum system. Among the various proposals 
for a physical realization of a qubit, specially designed 
superconducti ng c ircuits have been identified as promis- 
ing candidates .1321 Their potential is mostly founded on 
the high technological standards by which SQUIDS can 
be fabricated and controlled, as well as their promise re- 
garding array scalability .l^l Several types of Josephson- 
junction-based qubit designs, such as charge, phase, and 
flux qubits, have been proposed and explored in the lab- 
oratory.SI^lIzI Conditional gate operations have already 
been performed experimentally.'^'^ 

While preliminary results indeed look promising, con- 
siderable improvement in gate performance will be nec- 
essary to make these structures useful in larger arrays of 
quantum gates. Similar to all solid^tate-based qubits, 
Josephson qubits suffer from two main shortcomings: 
they are quantum two-level systems only within approx- 
imation and there is non-negligible coupling to the en- 
vironment. The former may result in leakage to non- 
computational basis states which affects the fidelity of 
quantum gates.^iSliniSIIIIIll The latter results in an un- 
wanted population relaxation and destruction of state 
superpositions, both being detrimental for quantum com- 
putation.LSMI 

In this work we study Josephson charge qubits and 
show how the conflicting requirements for suppression of 
both state leakage and decoherence can be accomplished 
within optimal control theory. In the flrst part, we out- 
line the model for the Josephson charge qubit includ- 
ing leakage and environmental interaction on which we 
base our study. Then, we formulate a cost functional 
for state-independent optimization of open quantum sys- 
tems based on the Kraus representation and process to- 



mography. This method can be readily applied to any 
open quantum system, including other qubit implemen- 
tations as well as multi-qubit gates and is not necessar- 
ily restricted to solid-state realizations. In the remain- 
der of the paper, we demonstrate the approach for the 
Hadamard gate implementation. 

II. JOSEPHSON CHARGE QUBIT 

The basic structu re of a Josephson charge qubit is 
shown in Fig. yj2l4Il2l The characteristic energy scales are 
the charging energy Eq of the superconducting island, 
the Josephson coupling energy Ej and the superconduct- 
ing energy gap A. If A is the largest energy, the problem 
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FIG. 1: (Color online) Superconducting charge qubit with 2 
Josephson junctions. The two junctions are realized by small 
insulating layers which divide a superconducting ring into a 
small island (red/gray) and a reservoir (thick black) .0 The 
qubit is driven by an external flux $ and a gate voltage V . 
The voltage is coupled capacitively (C) to the island. The 
Josephson junctions are characterized by a capacitance Cj 
and the Josephson coupling energy Ej. 

can be reduced to a situation where no quasiparticle ex- 
citation is found on the island and only Cooper pairs can 
tunnel through the Josephson junctions. Cooper-pair re- 
laxation on the island, as well as quasiparticles tunneling 
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from the reservoir onto the island, contribute, in princi- 
ple, to relaxation and decoherence. The latter takes place 
at time scales of the order (T|^)~^ ^ -^w where gr 
denotes the conductance of the Josephson junctions (in 
units of e^/h), 5^ is the quasiparticle level spaci ng in the 
reservoir and Nqp is the number of quasiparticlesl^^E^For 
typical charge-qubit architectures, T!^^ is of the order of 
10^"^ s to 10^^ s and, thus, on a time scale much larger 
than we are investigating in the present work. Cooper- 
pair relaxation within the island is even weaker and may 
be neglected, too. 

Within the charge basis { |?^)}, the Hamiltonian of the 
system, 



Hs - ^Ec [n - 7ic{V)Y 



may be written as, 



2E^j cos 



TT^jcose, (1) 



Hs = Y. {^^C [n - n,{y)f \n) (n| + 

n 

Ej{<^){\n + l){n\ + \n){n+\\)], (2) 

with Q (the phase difference across a Josephson junc- 
tion) being the conjugate variable to the number n of 
additional Cooper pairs on the island.l^ 



Ej{^) = -£;J}cos(7r$/$o) 



and nc{V) = VC/{2e) contain two independent physical 



control fields, V — Vn 



Vp and <&, where Vg and Vp 



denote a dc and a pulse gate voltage, respectively. <I> is 
the external magnetic flux. The latter can be used to 
tune Ej{^) between —Ej and O.I^SlBy proper choice of 
the dc gate voltage Vg, one can set the qubit "working 
point" to the charge degeneracy point, where the qubit 
is insensitive to charge fluctuations up to the first order. 
Using, next to the computational basis { |0) , |1)}, two 
adjacent leakage states | — 1) and |2), the effective "leaky 
qubit" Hamiltonian at the charge degeneracy point reads. 
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For our computation we consider a typical charge qubit 
with energi es Ec — 150 /leV and Ej = 35 fieV from the 
experiment .'^'i^'^Il xhe qubit can be controlled by tuning 
the external magnetic flux [<& = ^{t)] and the applied 
pulse gate voltage [Vp = Vp{t)]. 

In the derivation of the Hamiltonian [Eq. (ij], the 
loop-inductance of the superconducting ring has been ne- 
glected. For the parameters used in the present work, this 



approximation is justified if L <C ^q/ (^tt'^Ej) « 20 nH. 
We also assume that the Josephson coupling energies 
of the two junctions are identical. However, if E'j 
-Ej + (5-Ej for one of the junctions, Eq. (l3| turns into 



Ej{^) = -E0|(l + 7?)cos(7r|^) +Z77sin(^|^ 



with 



SE]'j 
2E" 



see e.g. Ref. (20( In order to avoid such an 

asymmetry, it is possible to substitute one of the Joseph- 
son junctions in Fig. [l] by a SQUID, which enables one 
to tune the Josephson energies to be equal.!^ 

In superconducting qubits the dominant dephasing 
mechanism (at low temperature) is attributed to noise 
which changes from l//to Ohmic behavior at frequen- 
cies of typically kBT/H)^ The microscopic origin of the 
1// noise is not fully understood yet but it is beHeved 
that it originates from background charge fluctuations. 
Ohmic contribu tions m ay result from intrinsic sources or 
from gate lines. '^^'^^'^ We map these fluctuations onto 
a bat h of harmonic oscillators coupling linearly to the 
qubit .1221211 The system is thus modeled by a Hamiltonian 
of the form. 

Hit) = Hs{t) + HB + Hi, HB = J2hL0kblbk, 

k 

Hi = hX®Y,9k [bl + bk) =X®T, (5) 



(3) where gk is the effective coupling constant of the spin- 



boson interaction and b\ and bk are the bosonic creation 
and annihilation operators for the mode with frequency 
LJk- The dynamics of system, bosonic bath, and the in- 
teraction between the two, respectively, is governed by 
the Hamiltonians Hs, Hb, and Hj, whereby the exter- 
nal control field e(t) to be optimized is contained in Hg 
only, thus assuming that there is no direct control over 
Hb and Hj. 

Apart from correlation functions, noise has to be clas- 
sified by its amplitude distribution (AD). For the spin- 
boson model with a thermal bath the AD is normal- 
distributed (Gaussian noise), which agrees quite well with 
experiment, especially at the short time scales with which 
we are concerned.'^^l Guided by the experimental results 
of Ref. [121 we construct the spectral density by choosing 
J{tu) = Ji/f{uj) + Jf{u}) consisting of a 1// and an Ohmic 
contribution,!^^ 

Ji//-«i///(4fcBT), Jf{Lu) = af/{2h)iu. (6) 

For thermal equilibrium, the noise spectrum Sr {ijj) asso- 
ciated with the operator F in Eq. (|5| in terms of J{ijj) is 
then (a tilde denoting the interaction picture). 



St (w) 



{r(t),f(i')} 



dr 



— 2hJ{ijj) coth 



2kBT 



{f(t),f(i')}e^™ 
T = t-t'. (7) 
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For hiv < 2kBT and Huj > 2kBT, Eqs. ^ and Q 
give the 1// noise spectrum, S'i/j(ti;) « ai/f/u, and 
the Ohmic spectrum, Si^{uj) « ajw, respectively. The 
strength of the charge fluctuations (proportional to 
ai/f), has been determined in experiment to saturate at 
~ (10~'^e)^ for temperatures lower than 200 mKl^^l The 
slope of the Ohmic contribution is given in Ref. 151 A 
plot of the specific noise spectrum which we use numeri- 
cally can be seen in the inset of Fig. [3^. 

Inspection of Eqs. Q reveals that simple control 
strategies for leakage suppression and minimization of 
decoherence are in conflict with one another. The most 
convenient way to minimize dephasing within the com- 
putational subspace would be a large absolute value of 
the Josephson coupling energy Ej{^) because it sets the 
minimum time (top) needed to perform unitary transfor- 
mations which incorporate rotations around the x-axis 
of the Bloch sphere (as needed for the Hadamard gate). 
The faster top, the lesser the effect of decoherence on 
the gate. However, in the present superconducting qubit 
architecture, Ej{^) is also responsible for coupling to 
the non-computational basis states | — 1) and |2). Due 
to the structure of X in Eq. ([s]), the net decoherence 
rate is enhanced when they participate in the dynamics. 
Thus, if decoherence-effects are present, the transfer of 
coherence to the leakage subspace is highly undesirable. 
In the unitary case, coherence can be transferred back to 
the computational subspace without loss (given sufficient 
control). These confficting requirements, as well as the 
complexity of an open quantum system, make the design 
of pulses which maximize the fidelity of the quantum gate 
a nontrivial task. 



III. STATE-INDEPENDENT OPTIMAL 
CONTROL 

In the context of quantum information processing, co- 
herent control of unitary operations within the qubit is 
of fundamental interest .'23 In contrast to the optimiza- 
tion of state-selective transitions, e.g. to maximize the 
probability of a certain pathway in a chemical reaction|2s] 
here we are interested in the optimization of the whole 
dynamical map. Hence, the qubit should perform a de- 
sired unitary operation, irrespective of its initial state. 
We call this task "state-independent" optimal control. 

We consider a Hilbert space which is a tensor product 
of the Hilbert space of the open quantum system (from 
now on denoted as system) TLs and the Hilbert space 
of the environment (which we will call bath) Tis, Ti. = 
T-Ls^T-Lb ^ {T-isi © Hs2) i^Ti-B- Within the system we 
distinguish between computational states on which the 
gate operation is specified, spanning Tisi, and "leakage 
levels" spanning Hs2-'^ 

The task of state-independent optimization can be 
stated as follows. The density matrix of the system, 
Ps{t), evolves in time according to the map, ps(0) i-^- 
Ps{t) = £f{ps(0)}, for which the superoperator £t is 



functionally dependent upon externally applied control 
fields e(t), i.e. £t = £t[f\- In order to maintain the pos- 
itivity of ps{t), the map £ has to be completely po sitive 
and thus can be represented by Kraus operators KmW^^ 



PS 



(t) = £,[e] {psm = KnAmpsmim), (s) 



with K„i[e]{t) depending on the propagator of the com- 
posite system, and hence on e{t). We now want to find a 
control field e*{t) for which, at some final time tj, £tf [e*] 
approaches the desired mapping £^1 as closely as possible. 
For quantum gate operations, for example, one would set 



£D{.} = UD{tf){.)UUtf). 



(9) 



where Uoitf) is the unitary operation to be executed 
within the gate. In the context of quantum informa- 
tion theory it is useful to formulate a cost functional 
within the language of process tomography (see e.g. 
Refs. 27l29l30l3ip to define a measure of how well such 
a, not necessarily unitary, desired operation has been ac- 
complished. We rewrite the mapping £ by expanding the 
Kraus operators, K^it) = J2n'^mnK„, with amn € C 
and Kn £ A, where A denotes a complete basis set of 
M X M matrices. 1221 M = Mc + Ml is the total number 
of orthonormal basis states, consisting of Mc compu- 
tational and Ml leakage levels. The final state of the 
quantum system, starting out in p{0), now reads, 



pitf) = £t^ {p{0)} = ^i^„,p(0)ifix™«(t/), 



(10) 



with Xmn — J2k '^kmCt'kn- ^'^^ example of the Joseph- 
son charge qubit described above, one has Ale — 2, and 
Ml = 2 with basis states { |0) , |1)} and { | ~ 1) , |2)}, 
respectively. Hence, we are dealing with a Af^ x Af^ = 
16 X 16 representation of x [see Figs.|2ja) and[2jd)]. The 
matrix x in Eq. ( 10 1 is usually termed process tomography 



matrix and in order to compute its elements, we choose a 
fixed set of operators {aj} = B (for simplicity we choose 
B — A) for which we determine the time-evolution with 
respect to the mapping £, i.e., 



^■j{if) =£tf Wj}. 



(11) 



For one-qubit operations with 2 leakage levels or two- 
qubit gates, the tensor product of Pauli matrices (g) Tj 
with i,j e {x,y,z,0} and tq = 1 is a convenient choice 
for A and B, which we use in present work. In ex- 



periment, Eq. (Ill corresponds to preparation in and 



subsequent time-evolution of suitable different initial 
states (~ aj) of the quantum system. Recently, process- 
tomograp hy me thods have been applied to superconduct- 
ing qubits. 1^2133] Because Eq. Q is a linear mapping, one 
can rewrite Eq. (11) by 



k 



(12) 
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The mapping Et^ is described completely in terms of co- 
efficients Cjk, whose experimental determination involves 
quantum state tomographyf23 The last step in order to 
calculate the process tomography matrix relates Cjk to 
Xmn- Therefore, we note that, combining Eqs. (10l-(12l, 



(13) 



To extract the coefficients Cjk we apply the scalar product 
(cTi, CTfc) = 6ik and get from Eqs. (|ll|-p^. 



Cji = (CTi,CT-,- (t/)) , 

= (gj, q-mgjQ-n) Xmn 



E 

rn.n 



(14) 

^irnjnXmn {if ) ■ 



Xmn{tf) can now be calculated by solving the system 
of linear equations {ai,aj{tf)) = J2m,n ^imjnXmnitf), 
which includes inversion of B. 
Defining the operator,!^ 



X ■ 



J2 iK ® K^)Xn 



(15) 



we formulate a simple cost functional, 

J^||Px-X^ir = tr{[Px-X^] [PX-X^]'}, 

(16) 

where P denotes the projector onto the Mc-dimensional 
computational Hilbert space and < J < Jmax — 2Mq. 
For the charge qubit example. 



P 



ij - ^ Si^kSj,k with i,j e {1,2, 

fc=6,7,10,ll 



16} 



The transformation given in Eq. ( 15 1 corresponds to 
stacking the columns of the density matrix ps from left 
to right on top of one another, i.e. ps col (ps), which 
yields a dimensional single-column vector, termed 
col (ps). For the present work, this transformation is of 
the following form. 



Pii 

P21 
P31 
P41 
Pl2 



PS 



Pll 


Pl2 


Pl3 


Pli 


P21 


P22 




P23 


P24 


P31 
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P33 


P3i 


_Pil 


PA2 


P43 


P44 



col {ps) 



P22 



P32 



Pi2 
Pl3 



P23 



P33 



Pi3 
Pli 
P24 
P34 
P44 



(17) 



where framed density matrix elements belong to the com- 
putational subspace. The process tomography matrix x 
changes accordingly, x ^ Xi given in Eq. (151. For 
further details see Ref. 28, J measures the norm dis- 
tance between the target operation x^ and the actual 
operation x executed at time t/ for control field e. For 
the Hadamard gate with 2 leakage levels we set. 



U*r,®UD, with Ud 



1 

7i 



0- 
110 
1-10 




A similar approach based on a super-operator formalism 
is given in Ref. jSH 



IV. QUBIT DYNAMICS 

Starting with the von-Neumann equation for the full 
system dynamics, ^p{t) = —j^ [H{t), p{t)], the time evo- 
lution of the operators crj(t) (with tilde denoting the in- 
teraction picture) is computed within a non-Markovian 
master equation in Born approximation,!^ 



dt 



(18) 



1 



dt'ivB { [Hi{t), \Hi{t'),a,{t') ® pb(0)]] }, 



using a numerical method described in Ref. |36] to deal 
with the kernel, which is non-local with respect to time, 
tre denotes the partial trace over all bosonic degrees of 
freedom. For ps(0) we assume a bath in thermal equilib- 
rium at T = 50 mK. We also introduce a sharp infrared 
and a continuous ultraviolet cutoff for the spectral den- 
sity, i.e., 

J{uo) ^ { Ji/f{Lo)e (w - A„) + Jfito)) e-^/^"". (19) 

We choose A^^ = 27r x 100 Hz and Auv = 2n x 
100 GHz (Ohmic noise in Josephson charge qubits has 
been measured up to frequencies of 100 GHz).E^ The 
lower bound for the infrared cutoff is determined by the 
data-acquisition time of the experiment. We have tested 
our optimal pulses with respect to different infra- and 
ultraviolet cutoff frequencies, see Figs. [Sj^a) and [Sj^b). 
When evaluating the trace of the double commutator in 
Eq. (18 1, one has to calculate bath correlation functions 



of the form 



f(i)f(i'))=trB{f(i)f(t')ps(0)}, 

which can be derived analytically using the spectral den- 
sity given in Eq. (19 I and 



V. RESULTS 

Our strategy is to use short control pulses of high am- 
plitude which move the qubit away from the degeneracy 
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point only briefly and are within current experimental 
capabilities. We employ a time-optimal control strat- 
egy which incorporates the final time of the gate oper- 
ation (tf) as an additional control parameter, compati- 
ble with experimental control field strength. The control 
fields are of the form, 



nc{t) 

m 

$0 



with |Ac/$,a;c/*,</3c/<i,,7c/$,t°/$| representing free pa- 
rameters. Additionally, we restrict {nc{t),^{t)} to start 
and end at the degeneracy point {uc = 1/2, $ = 0) 
by using an envelope function g(t) [dashed, gray line in 
Fig.[2][b)] and to satisfy other constraints imposed by the 
design of the superconducting circuit. There are other 
methods, such as driving the qubit by NMR-like tech- 
niques and adiabatic pulses, within which fidelities of 
~ 0.3 — 0.4 have been achieved.!^ Gate fidelities of ~ 0.4 
have also been reported for a CNOT gate implemented 
by two coupled flux qubits.-' Recently, single qubit op- 
erations with gate errors of 1 ~ 2% ha ve been realized 
within phase and transmon qubits.!^^!^ 



To optimize the cost functional Eq. ( 16 1 we use a par- 
allelized, constrained version of a differential evolution 
algorithm with 300 individuals per generation and 3000 
generations per optimization run.'^Hl The algorithm flnds 
an optimum for control flelds given in Figs.[2]jb) and[2|^c). 
We typically choose crossing probabilities (C) and scal- 
ing parameters {S) within a range of C £ [0.90, 0.96] and 
S G [0.75,0.85]. Variations in these parameters have a 
large impact on the convergence rate of the algorithm. 
All calculations have been performed on a 40-node SUN 
Linux Cluster. 

For t/ « 50 ps, gate fldelity (F) losses AF = 1 - F = 

{J[e*]/ Jmax)^^^ as low as « lO^'^ are achievable. A 
plot of the associated matrix-element deviations from 
the Hadamard operation, |(x — X^)^„|) is shown in 
Fig. [2|a). Losses from decoherence (mainly included in 
red/dark gray bars) dominate over those from leakage 
(mainly included in white bars) by about a factor of 15, 
largely due to Ohmic contributions in the noise spec- 
trum. To demonstrate the importance of including ad- 
ditional Cooper pair occupation on the superconducting 
island within the model, we optimized the gate without 
leakage and then used these control flelds to steer the 
qubit subjected to leakage. The value of the cost func- 
tional increases by about 4 orders of magnitude corre- 
sponding to an increase in AF by 2 orders of magnitude. 
Incorporation of leakage states | — 2) and |3) does not 
add signiflcant fldelity losses for optimal pulses given in 
Figs. [2ja) and[2jb). Therefore, expanding the leakage 
state-space beyond | — 1) and |2) is not necessary. 
We have also explored longer pulses {tf = 200 ps) and 



Short pulses ~ 50 ps 

(a) 

AbS[(A-A-°)mn] 



(b) 
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Longer pulses ~ 200 ps 

(d) 

Abs[<^-^°)„„] 



(e) 




50 100 150 200 

VPS 

FIG. 2: (Color online) (a,d) Deviation of the Hadamard pro- 
cess tomography matrix elements (x)m„ for optimized control 
fields [as given in Figs. 2]^b,e) and (c,f)] from desired values 
(y^) . Black areas denote matrix elements which are ir- 
relevant and, hence, arbitrary. They neither alter the qubit 
dynamics nor induce transitions out of the qubit subspace, 
i.e. leakage. For simplicity, they are plotted as zero [see also 
Eq. ( 17l]. Red (dark gray) columns show deviations of matrix 



elements which represent the gate operation within the qubit 
subspace, whereas white bars display leakage. (b,e) and (c,f) 
The strength of the optimal control fields lies within realistic 
values and the optimal solutions have a simple shape. 



do not start at the degeneracy point with respect to <& 
[the working point is now set at {ric = 1/2, <& — <&o/2), 
see Figs. |2jd)-|2jf)]. These pulses show an enhanced re- 
duction of leakage but they do lead to stronger deco- 
herence because of an increased pulse-duration. Also, 
moving the qubit 's working point away from the degen- 
eracy point causes higher sensitivity to fluctuations in 
the control flelds. Typical fldelity losses AF for 200 ps 

pulses are about 6 times higher than those for 50 ps ones 
(AF« 10^2)^ 

Time-optimal control inevitably leads to short pulses 
(typically of the order of 50 ps for present fleld intensi- 



6 



(a) 



J(AirAuv) 




2.0x10" 



1.6x10" 



1.2x10" 



8.0x10" 



4.0x10" 



(b) 
12 

11 

O 

^10 
< 

^ 9 



Air/S 



J(Ai,)atAu„=512GHz 



(c) 

16 

■V 12 
o 

S 8 

< 

^ 4 



J(Au„)atAir=1S 



gJ 2' 

Air/S"^ 



Au„/(10's"') 



FIG. 3: (Color online) (a) Variation in the cost func- 
tional with respect to changes in the infrared and ultravio- 
let cutoff frequency. We used the pulse sequences shown in 
Figs. [2|b) and ^c). Inset, solid line: Noise spectrum em- 
ployed for our calculations. Dashed and dotted-dashed lines 
represent 1// and / noise contributions, respectively, (b) For 
a given ultraviolet cutoff, the cost functional shows a log- 
arithmic dependence on Air, i.e. J (Air) ~ log(l/Air), as 
expected (see Ref. I23p . (c) For J{Auv) we obtain a logarith- 
mic dependence over a broad regime. However, when Ohmic 
contributions become relevant the cost functional increases 
significantly. 



ties) whose generation in experiment is a demanding task. 



However, picosecond electrical pulses can be produced 
e.g. by optoelectronic devices, such as photocunductive 
switches^ or by optical rectification of ultrashort optical 
pulses using nonlinear media (e.g. LiTaOs).'^ Ultrafast 
pulse-shaping methods have been discussed extensively 
in Ref. m 



VI. CONCLUSION 

By applying a newly developed, process-tomography- 
based optimal control theory for open quantum systems 
to a Josephson charge qubit, we have shown that one- 
qubit gates, such as the Hadamard gate, can be real- 
ized with remarkably high fidelity. The strategy has 
been to keep deviations of the control fields with respect 
to the degeneracy point as short as possible while per- 
forming the desired unitary operation in time-optimized 
fashion compatible with experimentally available con- 
trol field strength. A fully quantum mechanical descrip- 
tion based on experimental noise spectra has been em- 
ployed to model dephasing effects and additional non- 
computational basis states have been included to account 
for unwanted Cooper-pair occupation on the supercon- 
ducting island. Depending on the gate operation time, 
which has been treated as a variable, we could achieve fi- 
delities of the order ofF«l — 10~^tol — 10~^. In terms 
of the process tomography matrix this corresponds to er- 
rors of the order Abs [{x - X^)„„J ~ 10"^- 10"^. Thus 
we find that charge qubits can be made to perform at an 
equal level with current realizations of other Josephson 
qubits (such as phase, fiux or transmon qubits) which 
have been tested in experiment so far. 
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